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ABSTRACT 

Context. ZCMa is a binary composed of an embedded Herbig Be and an FU Ori class star separated by ~ 100 au. Observational 
evidence indicate a complex environment in which each star has a circumstellar disk and drives a jet, and the whole system is 
embedded in a large dusty envelope. 

Aims. We aim to probe the circumbinary environment of Z CMa in the inner 400 au in scattered light. 

Methods. We use high contrast imaging poiarimetry with VLT/NaCo at H and bands. 

Results. The central binary is resolved in both bands. The polarized images show three bright and complex structures: a common dust 
envelope, a sharp extended feature previously reported in direct light, and an intriguing bright clump located 0'.'3 south of the binary, 
which appears spatially connected to the sharp extended feature. 

Conclusions. We detect orbital motion when compared to previous observations, and report a new outburst driven by the Herbig star. 
Our observations reveal the complex inner environment of Z CMa with unprecedented detail and contrast. 
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1. Introduction 

Herbig stars are the earliest optically-revealed stage of medium 
mass star formation. They are usually found embedded in en¬ 
velopes rich in gas and dust, typically surrounded by circumstel¬ 
lar protoplanetary disks. Stellar multiplicity is an important ele¬ 
ment which has a strong impact on circumstellar evolution. Inter¬ 
estingly, the majority of stars are born in multiple systems ( [Kraus 


et al. 


2Q12| ). Z Canis Majoris (Z CMa) is an enigmatic binary sys- 


is surrounded by a huge infalling envelope ( [Beckwith & S argent j 
199 H[ Alonso-Albi et al.|2009| ), which can explain the v ery high 
accretion rate of the FU Ori star (~ 3 x 10“^ yr“^|van den 


Ancker et al.|2004[ ). Spectropolarimetric obseryations of ZCMa 


primary driyes a 3.6 pc jet ( [Poetzel et al.[[1989| ), and IFU ob¬ 
seryations with OSIRIS/Keck show that the secondary driyes a 
“microjet” (~ 400 au, Whelan et al.|2010 ). At a few hundred au 
south of Z CMa, ~ ' ‘ 

presence of a jet 
ing about 


Millan-Gabet & Monnier[ ( [2002[ ) reported the 


tern composed of an embedded Herbig Be star (p rimary) and an 
FU Ori (secondary) star separated by ~ 100 au ( [Koresko et al.[ 
1991[). Distance estimations to ZCMa range from 930 to 1150 
pc ( [Claria|1974[[Kaltcheya & Hilditch|2000|). The whole system 


reyealed the nature of the primary showing a spectrum rich in 
emission lines and largely polarized ( [Whitney et al.[[T99^ . A 
dusty enyelope or cocoon around the primary is most likely caus¬ 
ing this effect. The extinction (Ay ~ 10(5.2) for the primary 
(secondary), [Hinkley et al.|2013[ ) created by this cocoon can de¬ 
crease the brightness of the primary towards optical wayelengths 
(see also [Whitney & Hartmann|1993t[Szeifert et al.|2010[ ). The 


* Based on observations made with the VLT, program 094.C- 
0416(A). 


ike feature observed in scattered light extend- 
1000 au from Z CMa towards the southwest. It is not 
clear whether it is associated with a bipolar outflow, or a jet, or 
whether it is a completely independent structure. 

ZCMa’s light curve shows strong and moderate outbursts 
with variations both in duration and in strength. During quies¬ 
cence the secondary dominates the optical and near-infrared (IR) 
up to H band and the primary dominates at longer wavelengths. 
However, during the strong outburst episodes the Herbig star be¬ 
comes the dominant source at all wavelengths ( [van den Ancker[ 
[et al. [2004}[Hinkley et al.[2013] ). Between 2008 and 2009 Z CMa 
suffered a strong outburst (Amy ~ 2.5). While some observa¬ 
tional evidence favors an accretion burst from the primary as the 
source of this outburst ( [Benisty et al.|2010[ Hinkley et al.|2013[ ), 
other observations favor the formation of holes in the dust co¬ 
coon around the primary to explain the sudden increment in flux 
( Szeifert et al.[[201^ . Changes in the structure and/or optical 
thickness of the dust cocoon can naturally explain the outburst: 
the opening of new holes in the direction of the line of sight al¬ 
lows the light from the Herbig star to directly escape, resulting in 
a dramatic increment of the total flux of the system at optical and 
NIR wavelengths. [Canovas et al.[ ( |2012[ ) presented polarized dif- 
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Table 1. Observing log of the observations. Bad seeing frames do not 
contribute to the total exposure time listed below. All observations were 
taken during 19/01/2015. 


Target 

Band 

DIT 

[s] 

Tot. Exp Time 
[s] 

airmass 
[min, max] 

ZCMa 

H 

0.15 

957.9 

1.03, 1.15 



0.5 

166.0 

1.04,1.05 


Ks 

0.15 

603.0 

1.05 ,1.07 



0.5 

176.0 

1.07,1.10 

HD 52841 

H 

1 

160.0 

1.20,1.22 


Ks 

1 

160.0 

1.24, 1.26 


ferential images (PDI) at optical wavelengths, showing evidence 
of light escaping through a hole in the dust cocoon and cavities 
carved out by the two jets in the common dusty envelope sur¬ 
rounding the system. 

In this letter we present high contrast PDI observations of 
Z CMa at H and Ks bands with NaCo. Our results show the ex¬ 
cellent high contrast capabilities of this instrument after its im¬ 
plementation from VLT/UT4 to VLT/UTl. We find new polar¬ 
ized features down to an inner working angle of ~ 0'.'12, and 
resolve the two components of Z CMa in the two bands showing 
that the primary is driving a new outburst. Our results allow us 
to describe the complex inner environment of this young system 
with unprecedented detail. 


2, Observations 


We observed ZCMa and the comparison star HD 52841 with 
NaCo ( [Lenzen et al.|2003[|Rousset et al.|2003| ) on January 19th, 
2015, right after its commissioning at the VLT/UTl. The seeing 
was mostly stable (median ~ 0'.'87) during the observations with 
a few short episodes of bad (above I'.'IO) values and a median 
coherence time of tq = 3.5 ms. Z CMa was very bright {mh < 6) 
and the Adaptive Optics (AO) system delivered near diffraction 
limited images. 

Our observations were obtained with the NIR detector with 
the H and Ks filters. The S27 camera (27.05 + 0.1 maspx“^) was 
used in cube mode with the read-out set to Double_RdRstRd and 
the detector mode set to HighDynamic. We used individual ex¬ 
posure times of 0.15 s and 0.5 s for Z CMa, and 1 s for HD 52841. 
The frame loss rate in the short exposure images was minimized 
by windowing the read-out region to 512 x 520 px. The effect of 
bad pixels was reduced by dithering on the detector. HD 52841 
was all the time below the non-linear NaCo regime (10"^ counts). 
For Z CMa, the primary is saturated in all our images and sev¬ 
eral pixels surrounding both stars have values above the linear 
regime. We used NaCo in polarimetry mode combining a half¬ 
wave plate (HWP) with a Wollaston prism. In this setup the HWP 
rotates the polarization plane of the incoming light and the Wol¬ 
laston splits the light beam into two images with orthogonal po¬ 
larization states that are projected in different regions of the de¬ 
tector. A field mask avoids beam-overlapping. These two images 
are hereafter named as 4 (for the ordinary beam), and 4 (for the 
extra-ordinary beam). We took images with the HWP rotated by 
0°,-22.5°,-45° and -67.5° to reconstruct the Stokes parame¬ 
ters of the incoming light. The observations are listed in Table 


3. Data reduction 


We used our own pipeline outlined in |Canovas et al.| ( [2011[|2013| ) 
to process the observations. Here we briefly describe the major 
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aspects of the reduction. The images with the poorest AO correc¬ 
tion (obtained under seeing > I'.'l) were discarded in our anal¬ 
ysis. The rest of the images were dark subtracted, flat-field cor¬ 
rected and sky-background subtracted. Hot and dead pixels were 
flagged. The images were first aligned with a cross-correlation 
algorithm. We then refined the alignment up to an accuracy of 
0.05 px with a minimization method. In short, our algorithm 
finds the position that minimizes the standard deviation (cr) of 
the difference between the two images to be aligned using cubic 
interpolation. The polarized infor mation was extracted using the 
double-difference approach (e.g. [Hinkley et al.||2009t |Canovas| 
|et al.||2011| ). This way the 4 and 4 images are added and sub¬ 
tracted to produce an intensity and polarized image, respectively. 
When the HWP is at 0° the polarized image corresponds to the 
Stokes -i-Q image. Similarly, when the HWP is at 22.5°, 45° and 
67.5° the resulting polarized image corresponds to Stokes -fU, - 
Q, and -U, respectively. Instrumental polarization (IP) is usually 
corrected at this stage assuming that the central star is unpolar¬ 
ized (e.g., [Avenhaus et al.||MTi] ). Unfortunately, this is not the 
case for the two stars in the Z CMa system ( [Fischer et al.|1998| ). 
Additionally, the saturation of the primary and the extended de¬ 
tector area around the two stars with pixel values well above the 
non-linear regime of the NaCo detector complicate the polari- 
metric data reduction. We therefore measure the polarization de¬ 
gree in a 1-px thick (07027) contour around the non-linear region 
for each individual image. This region contains contributions 
from the IP and the interstellar polarization and reaches a max¬ 
imum of ~ 2.59^ in our observations. This quantity was sub¬ 
tracted from the individual images. To test for the effect of this 
correction we also reduced our images using different amounts 
of correction ranging from 1% to 4%. The non-linear pixels are 
flagged and not considered in our polarimetric analysis. The im¬ 
ages were median-combined to construct the final Stokes param¬ 
eters as 2 = (-rQ - (-2))/2 and U = (-\-U - (-U))/2. The (lin- 
early) polarized intensity ( Pi ) is described by Pj = A/2^~+~f4. 
The polarization angle indicating the vibration plane of the elec¬ 
tric field is Pe = ^arctan(UIQ). The total intensity is computed 
by adding-up the median-combined intensity images. Finally, the 
images were normalized to 1 s exposure time. 


4. Results 


The short-exposed images of Z CMa show the same structures 
as the long-exposure ones while having smaller inner working 
angle. Therefore we focus our analysis on the 0.15 s images. In 
what follows we adopt the most recent value of 930 pc for the 
distance to the source. 

The ZCMa binary is resolved in both H and Kg bands 
(Fig. upper panels). The primary component is saturated in 
both bands. To estimate the separation (p), position angle (PA), 
and difference in flux (AmH,Ks) of the Z CMa binary system is not 
straightforward because of the saturation of the primary and the 
small separation between both stars. We circumvent this prob¬ 
lem by fitting the primary’s point spread function (PSF) wings 
and the secondary’s full PSF. This way we derive a separation of 
p = O'.'l 12 + 07003 and a PA of 139.6° ±2.0°. The brightness dif¬ 
ference between the primary and secondary is noticeably larger 
at Kg than H band, with Amh = -LI +0.2 and Amks < -1-9, 
similar t o the values observed during an outburs t driven by the 
primary ( [Hinkley et al.|2009t Szeifert et al.|2010[ ). 


Comparing our results with those derived by [Millan-Gabet 
& Monni^ ( [2002[ ) using data from 2001 we detect an orbital 


Witzel et al. (20111 reports a maximum IP for NaCo at UT4 of 4%. 
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motion of APA = 10.6° +2.1° and Ap = O'.'OOS ± O'.'OOS in 14 
yr. Our results are in agreement with Millan-Gabet & Monnier| 
( 2002| ) who find a marginal change in p but a significant change 
of APA = 8?8 ± 1?5 when analyzing previous observations cov¬ 
ering a 11.2 yr baseline, deriving a APA rate of ~ 0?7 yr“^ for a 
circular orbit composed of a IM© and a 5Mq stars. 



Fig. 1. Top and bottom-left: Intensity images in logarithmic scale. AO 
artifacts are visible in the ZCMa images. HD 52841 is slightly elon¬ 
gated in the SE-NW direction. Bottom-right: The P/ image of HD 52841 
shows the noise that inevitably appears in the innermost O'.'l (indicated 
by the white circle) of the image as a result of imperfect speckle sup¬ 
pression. 


Our polarized images (Fig. show three complex features 
down to the inner ~ O'.'l2 from the stars in the two bands. First, 
we detect an approximately circular halo around the two stars 
(labeled as “1” in Fig.[^ left panel). The detector window trims 
the upper part of the image, where the halo is detected at sig¬ 
nal to noise (S/N) of ~ 8, suggesting that it extends beyond that 
limit. In the east-west direction the halo is detected up to ~ 0'.'8 
(~ 744 au) from Z CMa. Second, we detect a bright, irregularly 
shaped polarized clump at ~ 0'.'3 (~ 279 au projected distance) 
south from ZCMa (labeled as “2”). This clump is elongated in 
the east-west direction and it is separated by roughly O'.'l5 from 
the “kl” feature described by [Whelan et al.|p()10| ). Third, we 
re-detect with high S/N (~ 30) the s harp extended feature (la- 
beled as “3”) previously observed by [Millan-Gabet & Monnier] 
P002[ ) with Keck at J band (shown in white contours in Fig.|^ 
central panel). This feature extends down to 0'.'4 and appears 
spatially connected with the clump, and is likely to extend south 
beyond our field of view because of its high S/N (~28 at 1'.'5 
south). The relative polarized color (Fig.[^ shows that, in scat¬ 
tered light, feature “3” is remarkably redder when compared to 
the neutral (Amn-K, ~ 0) polarized halo or with the bright clump 
(AniH-Ks ~ ^•4)- Fig. 1^ also suggests that feature “3” is more 
extended in the east-west direction in Kg band than in H band, 
specially towards the south direction. In both bands the polar¬ 
ized flux along the feature decreases with projected distance to 
the star, but not uniformly. There is a local decrement in F/ at 
~ -1'.'3 south (see Fig. [2j central panel), also observed in in¬ 
tensity at J band in the Keck images. We find the same com¬ 


plex structures when processing the observations using different 
amounts of IP correction, therefore concluding that they are not 
artifacts. 

The polarization angle Pg at Ks band is shown in Fig. left 
panel. The vectors are plotted only in regions of the image with 
S/N >10 averaging over a 3-px width (~ 1.2 FWHM at Ks band) 
squared box. Overall the vectors follow an axi-symmetric pattern 
around the Z CMa binary system and they show the same ori¬ 
entation along feature “3”. The white arrows indicate polarized 
artifacts created by the AO system. In these regions Pg differs 
from the rest of the image, showing a non axy-symmetrical ori¬ 
entation. An extension towards the East seems to depart from the 
polarized clump (see Fig.|^ central panel). This feature is more 
likely an artifact, since it is located very close to one AO arti¬ 
fact, and it does not appear in any of the other images. A zoom 
around the stars is shown in the central and right panels of Fig. 

In these plots we only show Pg where the S /N > 30. 

Interestingly, there is a polarized filament connecting the 
clump (feature “2”) with the innermost regions of Z CMa. The 
bright pixels around the non-linear (set to 0) pixels are likely to 
have a contribution from the inherent noise observed in the in¬ 
ner O'.'l around the stars (see the Pj image of the comparison 
star at the bottom right panel in Fig.[T]). However, the two stars 
and their respective circumstellar dis ks (and cocoon arou nd the 
primary) are known to be polarized ( [Fischer et al.|1998} , and a 
large contribution from them and the polarized halo can explain 
the quasi-axisymmetric distribution of Pq, specially at Ks band. 


5. Discussion and Conciusions 


Our polarized observations show the complex and rich environ¬ 
ment of Z CMa. The nature of feature “3” remains unclear. It 
is not spatially close to the two jets detected by [Whelan et al.[ 
( 2010| ), and to date there are no detections of jets in continuum 
(while there are several detections in emission lines). Particularly 
intriguing is the apparent change in direction it follows before 
connecting to the bright clump, and the morphological similar¬ 
ities with broad cylindrical cavities recently observed in CO at 
late stages of star formation ( [Lee & Ho|2005[ [Arce et al.|2013[ ) 
and at lOyum around massive embedded stars ( [Li et al.[[2014| ). 
Although it is difficult to imagine how feature “3” can produce 
scattered light at < 1400 au (projected distance) from the stars, 
this feature could be explained by a large filament of dust with 
different scattering properties than its surrounding material. Al¬ 
ternatively, it could be created by light scattered ofiT the walls of a 
cavity created by the interaction of a molecular outflow with the 
surrounding material. We find no change in the shape of that fea¬ 
ture when comparing with the Keck images taken 14 years ago 
(Fig. central panel). Give n the youth of ZCMa (< 3 10^ yr, 
|van den Ancker et al.|2004| ) and the size of the polarized clump 
(feature “2”), a tantalizing explanation for this clump could be to 
relate it to an embedded, very young star. High resolution images 
at mid infrared wavelengths can help to disentangle the nature of 
this clump by providing evidence of a mid-IR excess (or lack 
thereof). 

We detect orbital motion on the binary system of ZCMa. 
Assuming a circular orbit composed of a IM© and a 5Mq stars 
and comparing with previous measurements, our data suggest an 
almost face-on orbit with a measurable orbital mo tion of APA = 
0?7 yr“\ in agreement with the values derived by Millan-Gabet 


& Monnier (2002). The derived relative photometry is consistent 
with a new outburst from the system. 

The 1 and 3mm continuum images obtained by [ Alonso-Albi 


et al. ( [2009 ) with IRAM-PdBI show a partially resolved structure 
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Fig. 2. Right and center: Pj images at H and Ks bands. Green contours indicate the position of the two stars. The two images are plot in the same 
scale, using square root scaling. The pixels falling in the n on-linear regime are set to zero. Left: the dotted and dashed lines indicate the position 
of large and micro jets, respectively ( [Whelan et al.|2010 i. Arrows “1”, “2” and “3” indicate the polarized halo, the bright clump, and the sharp 
extended feature, respectively. Center: Arro ws “a” and “d” shows an artifact a nd a local decrement in polarized flux, respectively. The white line 
contours the Keck J band image presented in jMillan-Gabet & Monnier| ( |2002j ). Right: Relative polarized color computed as [H] - [KJ. The white 
cross indicates the brightest point of feature “2” at Ks band. 
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Fig. 3. Left: Pq over polarized intensity at K^-band. The length of the vectors is arbitrary. The image is plot in lineal scale, stretching the scale to 
highlight faint features. AO artifacts are indicated by the arrows. 


elongated in the N-S direction that could be related to the obser¬ 
vations here presented. In any case, images at (sub-) mm wave¬ 
lengths with much higher spatial resolution are urgently needed 
to study the cold-dust counterpart of our observations to com¬ 
plete our picture of Z CMa and to understand the earliest evo¬ 
lutionary stages of the binary systems. Finally, our results show 
that NaCo’s imaging polarimetric capabilities are at its best, al¬ 
lowing us to probe the innermost regions of Z CMa 
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